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Abstract 



There is a wealth of applied problems that can be posed as a dynamical system defined on a network 
with both attractive and repulsive interactions. Some examples include: understanding synchronization 
properties of nonlinear oscillator;, the behavior of groups, or cliques, in social networks; the study of 
optimal convergence for consensus algorithm; and many other examples. Frequently the problems involve 
computing the index of a matrix, i.e. the number of positive and negative eigenvalues, and the dimension 
of the kernel. In this paper we consider one of the most common examples, where the matrix takes the 
form of a signed graph Laplacian. We show that the there are topological constraints on the index of 
the Laplacian matrix related to the dimension of a certain homology group. In certain situations, when 
the homology group is trivial, the index of the operator is rigid and is determined only by the topology 
of the network and is independent of the strengths of the interactions. In general these constraints give 
upper and lower bounds on the number of positive and negative eigenvalues, with the dimension of 
the homology group counting the number of eigenvalue crossings. The homology group also gives a 
natural decomposition of the dynamics into "fixed" degrees of freedom, whose index does not depend on 
the edge-weights, and an orthogonal set of "free" degrees of freedom, whose index changes as the edge 
weights change. We also present some numerical studies of this problem for large random matrices. 

1 Introduction 

1.1 Problem Formulation 

There are many applied problems that can ultimately be reduced to the question of understanding a dynam- 
ical problem on a network or graph. In these applications it is often important to understand the dynamical 
behavior of the evolution in terms of the topological properties of the graph. 

In this paper we consider a simple, connected, undirected edge-weighted graph T = (V(T), E(T)) with 
vertex set V^r) and edge set E(T). For each edge in E(T) connecting vertex i with vertex j we associate a 
weight jij , which is assumed to be non-zero but may take either sign. If there is no edge connecting vertices 
i and j the weight 7^ is understood to be zero. 

For such a graph we define the signed Laplacian matrix £(T) by 



Note that £(T) is symmetric, so all eigenvalues of £(T) are real. If the weights are all positive, 7^ > 0, this is 
a standard graph Laplacian and as such the spectrum of £ is well-understood: £ is a negative semi-definite 
matrix, with the dimension of the kernel equal to the number of connected components of the graph T. 
In many applications the weights are not guaranteed to be positive. In this case the matrix is no longer 




(1.1) 



definite, and we are interested in determining the number of positive, zero, and negative eigenvalues of 
£(r), which we denote as n + (T), n (T), n_(F) respectively. 

In particular, we would like to determine topological bounds on these indices, i.e. conditions that depend 
only on the arrangement of signs of the edge weights 7^, and not on their magnitudes. We are aware of at 
least five applied problems that motivate this question: 

1 . Stability of Fixed Points Consider any dynamical system defined on the network. Specifically given a 
graph T, and symmetric coupling functions 92^ ( • ) = fji(-) associated to edges in the graph we define 
a dynamical system as follows 

4rXi = Fi(x) := uj, + ^ <Pij(xj ~ x % ). (1.2) 
jer 

A well-studied example of this type of dynamical system is the Kuramoto oscillator [2 29, 30. 38)|40), 
where we choose (fij(-) — jij sin(-). 

To compute the stability index of a fixed point for the system < |1.2) , i.e. a vector x with Fi (x) = for 
all i, we need to determine the index of the Jacobian J, where 



J2kfik( x k-x l ), i=j. 



The Jacobian J is a graph Laplacian of the form thus, determining the stability indices for fixed 
points of ( |1.2) is equivalent to the problem studied here ||4) [3T[|33}|42)[43) . When studying this dy- 
namical system, the first object of study is always the stable points. But, for example, if we consider 
such a system perturbed by a small stochastic process, then gaining a qualitative understanding of 
the dynamics requires that we identify all of the 1-saddles, i.e. those points which are unstable but 
with exactly one unstable direction (e.g. see [11 1). 

The (p' i: j terms in J can be of either sign, making it natural to consider the case of a graph Laplacian 
with arbitrary signs on the weights. For a generic choice of coupling functions, ip'^ (xi — Xj ) is non-zero 
for all i, j, implying that the graph determining J and the graph defined by the original interactions 
in \l.2\ have the same underlying topology. 

2. Stability of Neural Networks If we consider any neural network system with both "positive" and 
"negative feedbacks", then the stability analysis reduces to a eigenvalue problem similar to < |1.1) , In 
biological applications the interaction strengths are both difficult to measure experimentally and (due 
to neural plasticity) very changeable. This makes it impractical to estimate the magnitudes of 7^ with 
any degree of reliability. On the other hand the nature of the interaction (excitatory or inhibitory) is 
anatomical and, in general, will not change. Thus we have reliable experimental data on the signs of 
various connections, but very little reliable information on the magnitudes. This makes the idea of 
understanding the extent to which the dynamics is determined by the topology of the network a very 
attractive one. 

3. Convergence of Consensus Algorithms It was pointed out by Xiao and Boyd [44] that for some 
graphs, the optimal choice of weights for the convergence of a consensus algorithm uses negative 
weights. In particular, given a weighted graph V = {7^}, define the linear discrete-time consensus 
algorithm given by 

z<(t + l) := T^jijXjit). (1.3) 



The question posed in [44] is how one might optimize the choice of 7^, given the underlying graph 



topology and the constraint that 1 be a stable fixed point of 1 1.3 1, to obtain the most rapid convergence 



to consensus (to wit, to make the Lyapunov exponent of 1 1.3 1 as small as possible). It was observed 
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there that there exist examples of graphs where the optimal choice involves negative weights — that is 
to say, that there are choices of weights that make | |1.3) converge more rapidly than the most rapidly- 
mixing Markov chain. This observation inspired a deluge of work (examples include [17 22 .35 ,45 1) on 
this fast convergence problem and it has been observed that the need for negative weights is typical 
in many contexts. 

Of course, one obvious constraint on determining the optimal choice of weights is that n + (T) = so 
that we obtain convergence at all. In the work mentioned above, this was always obtained by solving 
a semidefinite programming problem over the set of all weights associated to a particular unweighted 
graph T. The results in the current work give topological (i.e. weight-independent) bounds for n+(F). 

4. Clustering in Social Networks Questions of this sort arise in the study of social networks. In the 
classic work of Hage and Harary [16], matrices of the form modeled the interaction of tribal 
groups within an alliance in New Guineaj^] In this work the underlying graph has sixteen vertices, 
representing the different tribal units, with edges represent relations between different tribal groups. 
These relations can be friendly ("rova") or antagonistic ("hina") corresponding to 7^ > and 7^ < 
respectively. Another anthropological example is the SlashDot Zoo, a social network associated to 
the website SlashDot |1|. On this website participants can label each other as friend or foej^] In the 
context of social network models the index n + (T) indicates the tendency of the network to separate 
into mutually antagonistic subgroups. 

5. Resistor Networks with Negative Resistance One intuitively appealing way to think about a signed 
graph Laplacian is as a network of resistors where some of the edges have negative resistance. While 
perfect negative resistors cannot be implemented with passive components there are many nonlinear 
components whose current-voltage curves are non-monotone, and have a region where the current 
is a decreasing function of the voltage. In the region where the current is a decreasing function of 
current the component has negative (differential) resistance. Some examples of such components 
include tunnel and Gunn diodes, neon lamps, and certain kinds of tubes. Circuits built from these 
components often exhibit multi-stability, and the stability of a given solution is governed by a system 
of the form Alternatively a negative resistance can be implemented using active circuit elements 
with gain, such as an operational amplifier (op-amp). 

On a somewhat different note we should remark that signed graphs and signed Laplacians arise very 
naturally in knot theory. One classical result of this kind is a procedure for associating a quadratic form to 

14) . In Goeritz's construction one considers a planar projection of the 
knot, with regions of plane alternately colored black and white. Two regions of given color are connected 
by an edge if they share a crossing. The weight attached to the edge is +1 if the crossing is left-handed and 
— 1 if the crossing is right-handed. The reduced determinant — the product over the non-zero eigenvalues 
of the associated graph Laplacian — can be shown to be a knot invariant. Related results include a con- 
struction by Kauff man of a Tutte polynomial for signed graphs 1 25 ] which specializes to known invariants 
like the Jones [21 1 and Kauffman bracket polynomials |24|. This construction was further generalized to 
matroids by Zaslavsky [49]. Signed graphs and their Laplacians have also been studied in the graph the- 
ory community independent of their connection to knot theory (see [46-48 1 and following, also |18|) but 
the the reader should be aware that there are several different generalizations of the Laplacian to the case 
where negatively weighted edges are allowed. In one variation definition the diagonal entries are taken 
to be minus the sum of the absolute values of the edge weights. In this case the matrix is typically not 
zero-sum, and is negative semi-definite. Clearly the spectral questions are different and there is no obvious 
correspondence between the two. 

: It should be noted that Hage and Harary were concerned with the question of balance in signed graphs, a very different question 
from the ones we consider here. 

2 In this example the edges are directed, leading to a non-symmetric Laplacian matrix. As we do not consider non-symmetric 
graphs in this work, we will need to consider a symmetrized version of the SlashDot Zoo in the applications below. However, many 
social networks (e.g. Facebook) are symmetric social networks by definition. 



a knot originally due to Goeritz 1 13 
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To make clear the kinds of questions we want to ask, and the kind of phenomena we would like to 
understand, we begin with an example of two networks. 




o 




/3 
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Figure 1: Two graphs on 15 nodes. Here we are represented positively weighted edges with solid black 
lines, and negatively weighted edges by dashed red lines. 

Example 1.1 Figure [I] depicts two signed graphs representing networks, in this figure dashed lines (red online) 
represent negative/excitatory fhina interactions and the solid line represent positive/inhibitory frova interactions. Each 
of these graphs has fifteen vertices, fifteen hina edges and fifteen rova edges. The dynamics, however, are typically very 
different. We claim that there is a marked difference between the way mathematicians tend to view these problems and 
the way other network scientists view them, and that these differences are summarized in the graphs above. The first 
graph is more symmetric, having as its automorphism group (in the absence of weights) the dihedral group Di 5 . The 
automorphism group of the second graph is trivial, consisting of only the identity. Since symmetry reduction is one 
of the most common and most powerful techniques in applied mathematics this suggests that the dynamics of the first 
network should be easier to understand. Many network scientists, however, would claim that the second network is 
simpler and easier to understand. They would argue that this network can clearly be separated into four functional 
units. The largest functional unit consists of vertices one through seven, which interact through mutual inhibition. 
Hanging off of this structure are three small units: the first consisting of vertices one and eight through twelve; the 
second of vertices two, thirteen and fourteen; and the third of vertices three and fifteen. These vertices in these smaller 
units act on each other via excitation and couple to the main unit through the shared vertices. The goal of this paper is 
to mathematize this intuition about the structure of the network. We will show that, at least for the kinds of questions 
asked in this paper, the second point of view is the more fruitful and amounts to an observation about a particular 
homology group of the graph. The first network has the property that a certain homology group has the maximum 
possible dimension. This, as we will show, implies that the spectrum of the Laplacian is essentially arbitrary: it has 
one zero eigenvalue, and the signs of the other fourteen eigenvalues can be chosen arbitrarily with an appropriate 
choice of weights. In the second network, on the other hand, the analogous homology group is trivial. This implies 
that the spectrum is rigid: the Laplacian always has eight negative eigenvalues, one zero eigenvalue, and six positive 
eigenvalues regardless of the choice of weights on the edges¥\ 




In this paper we give the best possible bounds on n_(T), Uo(r), n+(T) involving only topological in- 
formation — connectivity of the graph and the sign information on the edge weights. For a graph with 
N, vertices the difference between the upper and lower bounds is an integer that can vary between and 

3 Note that there is no obvious relationship between the size of the automorphism group and the size of the homology group. In 
these examples they are both large or both small, but it is easy to construct examples where one is small and the other large. 
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N — 1, depending on the topology of the graph. This integer represents the dimension of a certain homol- 
ogy group, and counts the number of possible eigenvalue crossings from the left to the right half-plane. 
The examples above represent extreme cases. In the first example the homology group has fourteen gen- 
erators and there are fourteen possible eigenvalue crossings. In the second example the homology group 
is trivial, consisting only of the zero element, and there are no eigenvalue crossings. This homological con- 
struction gives a natural splitting of the vector space into a "fixed" subspace, where there cannot be an 
eigenvalue crossing, and a "free" subspace, where all of the eigenvalue crossings occur. We also show that 
these bounds are strictly better than those implied by the Gershgorin theorem. Finally we will conclude 
with some numerical experiments and examples. 



2 Main Theorem 



2.1 Preliminaries 



The main theorem, Theorem 2.8 gives tight upper and lower bounds on the number of positive, negative, 
and zero eigenvalues. To state the main theorem of the paper, we first present a few definitions. 

Definition 2.1 Given a graph V, we define the two subgraphs T + (resp. r_ ) to be the subgraphs with the same vertex 
set as T (V(T) = V{T + ) = V(T-)) together with the edges of positive (resp. negative) weights: 

E(T+) = {ee E(T)\ 7ij > 0}, E(T_) = {e e E(T)\ 7ij < 0}. 

Moreover, these graphs inherit the weighting from the original T, i.e. 

( r +)y = max (7y > 0), (r_) y = min(7 u -, 0). 

We further adopt the following definitions and notations: 

• We use the notation c(-) to represent the number of components of a graph. 

• We let , i e (1 . . . c(T + )) denote the i th component ofT + , and similarly T_. 

• Given a graph V we let STiT) denote the set of all spanning trees of G. More generally we let STfc(r) denote 
the set of all spanning trees ofT having exactly k edges in T_. Note that U^JT^STk = ST and STkHSTk' = 
ifk £ k'. 

• We define the flexibility of a weighted graph as the number 

r(T) :=mr)|-c(r_)- c (r + ) + i. 

Ifr(T) = 0, then we say that T is rigid. We show below that the flexibility is always a non-negative number. 

• As mentioned above, for any weighted graph T, we define the three indices no(r), n_(F), n + (T) as the number 
of zero, negative, and positive eigenvalues of £(T). We will occasionally use the notation ind(r) to represent 
the triple ind(T) = (ra_(r),n (r),n + (r)) . 

Remark 2.2 Any graph Laplacian as defined in has all row sums equal to zero, so that C(T)1 = 0, and one 
necessarily has no(T) > 1. It is well-known [ 7, 1 0[ that if all the weights 7^ > 0, then the graph Laplacian is negative 
semi-definite, with n (T) = c(T) and thus ri- = \ V(T)\ — c(T). In particular, ifT is connected with positive weights, 
then 

n+(r)=0, no(r) = l, n_(r) = \V{Y)\ - 1. 

This is no longer true when the edge-weights are allowed to be negative — the Laplacian matrix of a connected graph 
can have multiple zero and positive eigenvalues. 
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Lemma 2.3 Every signed graph T satisfies the inequality 



e(r + ) + c(r_)<|v(r)| + c (r). (2.1) 

IfT is connected, writing \ V(T)\ = N gives 

c (r+) + c(r_) < n + 1. 

From f/zz's, it follows that the flexibility t(T) of any graph is a non-negative integer. 

Proof. We first note that if T_ contains a subforest T with t edges, then c(r_) < \y\ — I. 

We now define T = (V , E) to a quotient graph of Y, where V are the connected components of T + and 
(a, b) e E iff there is at least one edge in T_ from a vertex in component a to a vertex in component b. 

Since T is connected, so is T. Consider any spanning tree T of T. By definition, this contains c(T + ) — 1 
edges, since it is a tree on c(r + ) vertices. Now consider this tree "lifted" into T where, for every edge in T of 
the form a f-> b, choose one edge in E 1 that connects component a to component b. This must be a subforest 
of r_, since it cannot contain any cycles. Therefore c(r_) < \V\ — (c(r + ) — 1) and we are done. 

Finally, if T is not connected, let r' 1 ) , ... , r( c ( r " be the connected components of T, and define 1+ , r_ 
in the obvious manner. By assumption, we have that 



c(r<?) + c (rw) < 



p(i) 



It is not hard to see that 



and from this (12. 1 b follows. 



C (r) c (r) 

^ C (r«)^ c (r + ), |r| = ^|rw 

i=l i=l 



Remark 2.4 Since t(T) is a counting number, the obvious next thing to determine is what it is that it counts. We 
will show in Section 2.3 from the consideration of a particular Mayer -Vietoris sequence that t(T) is the dimension of a 
certain homology group. This will be useful because it gives a relationship between the flexibility of a graph and cycles 
of a certain type. However, it is convenient to have the above proof, which is self-contained and purely graph-theoretic, 
at the current time. 



The main machinery that we will use in this paper is the celebrated Kirchhoff matrix tree theorem [6 26 
[27) . To state it, we first present some notation, which essentially follows that of Tutte pl| : 

Definition 2.5 • Let The a tree, then we define n(T) to be the product over the edge weights in the tree 

<T)= [] la- (2-2) 



• Let r be a weighted graph with \V(T) \ — N, and C(T) be its graph Laplacian. We know that C(T) has a zero 
eigenvalue, and therefore det(£(r)) = 0. Order the n eigenvalues of L{T) so that Xi = 0, then we define 

(_-i)N-l N 

• A/f ( r )= N Ii X * ( 23 ) 

i=2 

In other words M(T) is (up to the multiplicative prefactor) the linear term in the characteristic polynomial of 
the Laplace matrix. Note that A4(T) ^ OiffO is a simple eigenvalue of C{T). More generally, the multiplicity 
of a zero of M{T) is one fewer than the multiplicity of zero in the characteristic polynomial of C(T). M(T) is a 
reduced determinant which is zero iff C(T) has a non-simple zero eigenvalue. 
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With this notation the Kirchhoff matrix tree theorem can be stated as follows: 

Lemma 2.6 (Weighted Matrix Tree Theorem) Let T be a connected, weighted graph, and ST(T) the set of all 

spanning trees ofT. Then 

M(T)= < T )- (2-4) 

TeST(r) 

Remark 2.7 This is Theorem VI.29 in the text ofTutte pl| : a proof is provided there. Notice that if all of the edge 
weights are non-negative, then the sum in \2A\ is a sum of positive terms. This is an alternate proof that the kernel 
of a graph Laplacian with positive weights is simple for a connected graph. However, once we allow negative weights, 
the sum on the right-hand side can have cancellations, and this is the major difficulty in understanding the spectral 
properties of graphs with negative weights. 

2.2 Statement and proof of main theorem 

We begin by stating one of the main results of this paper. 

Theorem 2.8 Let T be a connected signed graph, and n_(T), n (T), n + (T) be the number of negative, zero, and 
positive eigenvalues respectively. Then for any choice of weights one has the following inequalities: 

c(r+)-i <n+(r) <7V- C (r_), 

c (r_)-i<n_(r)<iv- c (r + ), (2.5) 

1 <n (r) < N + 2 - c(T_) - c(T+). 

Further these bounds are tight: for any given graph there exist open sets of weights giving maximal number of 
negative eigenvalues 

n+(T) = c(T+) - 1, 

n_(r) = 7v- c (r+), 

no(T) = 1. 

as well as open sets of weights giving the maximal number of positive eigenvalues 

n+(T) = N - c(T_), 
n_(r) = c(T_)-l, 
no(T) = 1. 

Remarks 2.9 Notice that in each inequality in \2.5) the difference between the upper and lower bound is exactly 
the flexibility of the graph t(T) . This shows that t(T) counts the number of potential eigenvalue crossings. More 
precisely the theorem shows that there are c(T + ) — 1 eigenvalues which are always negative, c(F_) — 1 eigenvalues 
which are always positive, and r(F) = N + 1 — c(T_) — c(T + ) eigenvalues depend on the choice of weights. For rigid 
graphs there are no eigenvalue crossings and the index is fixed regardless of the choice of weights (thus inspiring the 
terminology "rigid"). 

In general the index will be constant with n (T) = 1 on open sets, separated by codimension one sets where 
n (T) = 2. One expects that n (T) > 2 only on sets of higher codimension. 

In the context of network models such as \1.2\ described above, the theorem shows that a necessary condition for 
stability of a phase-locked state is that between any two oscillators there exists a path such <p'ij{xi — Xj) > Ofor all 
edges on the path. While the necessity of such a condition seems physically obvious we are unaware of a previous proof 
of this. 

The main idea of the proof is to consider a one-parameter family of weighted Laplacian matrices. We 
then compute a polynomial associated to the graph, the zeroes of which detects eigenvalue crossings. We 
then show that this polynomial has exactly r roots in the positive half -line, and that the multiplicity of these 
roots is equal to the number of eigenvalues crossing from the left to the right half -plane at that parameter 
value 
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Definition 2.10 Given a weighted graph T we define a one-parameter family of weighted graphs T(t) as follows: the 
weights ofT(t) are related to those ofT by 



Jij(t) : 

or, more compactly, T{t) = T + + £T_. Obviously T = T(l 




We also recall Lemma 2.6 the weighted matrix tree theorem, and define 



M{T{t))= < T )- @ 

Tesr(r) 



The next observation is that A4(T(t)) is a polynomial of a very special form: 
Lemma 2.11 M(T(t)) is a polynomial in t which takes the following form 



jv- c (r_) 
fc= c (r + )-i 



where the coefficients a& are given by 



Tesr k (r) 



All of the a k appearing in \2.6\ are nonnegative; moreover, the first and last coefficients, a c ^ r+ - ) _i and ajv_ c (r_)/ 
strictly positive. 

Proof. Clearly Ai(T(t)) is a polynomial in t, since it is given by sums and products of terms each of 
which is at most linear in t. Since there is one power of t associated with each negatively weighted edge 
and ST = U%~£ST k it follows that 

JV-l 

M(T(t)) = ^-t) k - 

fc=0 

with defined as above. 

Next we note that STk is empty for k < c(T + ) — 1 and non-empty for k = c(T + ) — 1. To see this note that, 
since T + has c(T + ) components we need at least c(T + ) — 1 negative edges to connect them, and that there is 
at least one way to construct a spanning tree with c(T + ) — 1 negative edges: we first construct a spanning 
tree on each of the c(T+) — 1 components of T + using positive edges, and then connect them with c(T + ) — 1 
negative edges. This can always be done since the graph is assumed to be connected. Since a c (r + )-i is a 
sum over a non-empty set of positive terms it is positive. The upper bounds follow from the dual argument: 
reversing the roles of T + and T_ shows that any spanning tree must have at least c(T_) — 1 positive edges. 
Since any spanning tree has exactly N —1 edges there are at most N —I — (c(T_ ) — 1) = TV — c(T_ ) negative 
edges. 



Remark 2.12 We note that the polynomial M.(T(t)) is strongly reminiscent of other graph polynomials such as 
the chromatic, rank and Tutte polynomials, which have definitions in terms of sums over spanning trees. We will 
show later in the paper that M{T{t)) satisfies a contraction-deletion relation similar to that satisfied by other graph 
polynomials. 

Lemma 2.13 The roots of the polynomial M{T(t)) are real and non-negative. 
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Proof. This follows from the observation that the roots of the polynomial are exactly the eigenvalues of 
a generalized symmetric eigenvalue problem. Note that a root of the polynomial Ai (T(t)) corresponds to a 
solution of 

C + v = —tC_v 

where v can be assumed to be orthogonal to(l,l,l,...,l).A standard result in the theory of the generalized 
symmetric eigenvalue problem (gsep) is that a sufficient condition for the problem to have all real eigenval- 
ues is that there exists a linear combination of L + and L_ that is strictly positive definite. For completeness 
we give a short proof of this here. First note that the gsep Av = XBv has real eigenvalues if either A or B 
is strictly positive definite. If B is strictly positive definite then the above problem is self-adjoint under the 
inner product (v, v) — v l Bv, proving reality of the eigenvalues. If A is strictly positive definite a similar 
calculation holds with A replaced by A . Next note that there is an equivariant action of GL(2, R) on the 
generalized symmetric eigenvalue problem, as follows: if (A, B) is a pair with eigenvalue A Av = XBv then 
for any scalars a, /3, 7, 6 with aS — f3j 7^ the pair (aA + f3B, jA + SB) has an eigenvalue /i = and the 
same eigenvector v 

(aA + j3B)v = [ifrA + SB)v. (2.7) 

Thus if there exists a positive linear combination of A and B then the gsep Av = XBv has only real 
eigenvalues. In the case of L + and L_ it is clear that (for instance) L + — L_ is strictly positive definite 
on (1,1,1,...,1), as it is the graph Laplacian for a connected graph with positive weights. Therefore 
C + v = —tC-v has only real roots. 

For more information on the gsep see the review paper of Parlett [36], Theorem 1 in the paper of Craw- 
ford 18 9 1 or chapter IX §3 of the text of Greub |15|, for a proof (due to Milnor) of a somewhat stronger 
result. 

That none of the eigenvalues can be negative is clear from the form of the polynomial: since the coeffi- 
cients alternate in sign p(t) is obviously non-zero for t < 0. 

■ 

Remark 2.14 The fact that Pr(t) has only real roots implies, via Newton's inequality, that the sequence {ak}^ =0 is 
log-concave 

^ 2 

flfe+ia/c-i S ak- 
in the special case where the weights of the edges are all ±1 then the coefficients au are integers which count the 
number of spanning trees ofT having exactly k edges in r_. A large number of other combinatorial sequences share 
this property. See the review papers of Stanley [39] or Brenti [3] for details. The analogous problem of the log-concavity 
of the coefficients of the chromatic polynomial, a much more difficult problem, was a long-standing conjecture that has 
recently been established by Huh 119 20|. 

Next we show that the multiplicity of the zeroes of the polynomial A4(T(t)) is equal to the dimension of 
the kernel of 1,1,. First a preliminary lemma: 

Lemma 2.15 The eigenvalues Xi(t) corresponding to eigenvectors orthogonal to (1, 1, 1, . . . , 1) are non-decreasing 
functions oft that cross zero transversely: ifXi{t) = then A-(i) > 0. 

Proof. The fact that Aj(t) are non-decreasing follows immediately from the fact that £ + is a positive 
semi-definite matrix. The transversality follows from a perturbation argument and elementary topological 
considerations as follows. If A^ (t) vanishes at t = t* then from degenerate perturbation theory |23| we have 
that A^(i*) is equal to one of the eigenvalues of the matrix 

£+lkcr(£(i*)) 

The matrix £ + is positive semi-definite, and thus A^(t*) > 0. To see the strict inequality we note that in 
order for £+|ker(.c(t*)) to have a zero eigenvalue there is necessarily a vector in ker(£ + ) D ker(£_). Such a 
vector would be constant on components of T + and T_ and thus, by the connectedness assumption, on 
all of T. The only vectors in ker(£ + ) n ker(£_) are thus multiples of (1, 1, 1, ... , 1), and so any other zero 
eigenvalue must cross through the origin transversely. ■ 
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Lemma 2.16 The dimension of the kernel of C(t*) restricted to (1, 1,1,..., I) 1 - is equal to the multiplicity oft* as a 
rootofM{T(t)). 

Proof. The polynomial M(T(t)) = Xi(t) where, from the above, the Xi(t) have only simple roots. 

Thus the multiplicity of a root of A4(T(t)) is equal to the number of Aj(t) that vanish there. ■ 



We are now in a position to prove Theorem 2.8 



Proof of Theorem 2.8 The matrix £(T_) is positive semi-definite so that the eigenvalues of C(T(t)) are 



non-decreasing functions of t. For t = the matrix C(t) is a graph Laplacian with c(T+) components, so it 



has a c(T+) dimensional kernel and an TV — c(T+) dimensional negative definite subspace. By lemma 2.15 



c(r + ) — 1 of these zero eigenvalues cross transversely into the positive half -line, so for t small and positive 



the index of C(T(t)) is (N— c(T+), 1, c(T+) — 1). By lemma 2.13 the crossing polynomial has exactly r roots on 
the open positive half -line, and each root of the crossing polynomial corresponds to an eigenvalue crossing 
from the left half -line to the right, so one has exactly r eigenvalue crossings. This gives 

c(r+) - i < n+ (r) < n — c (r_) 
c(r_) - 1 < n_(r) < n - c(r+). 

When t is small and positive the lower holds for n + (T) and the upper bound for n_(T), and vice-versa 
when t is large. ■ 
In the next section we consider more detailed asymptotics on the eigenvalues in the limits ( < 1 and 

<> 1. 



2.3 Topological characterization of the flexibility 

The quantity r is a measure of the flexibility or rigidity of the network dynamics: it is a measure of the 
number of eigenvalues that can cross from the right half-plane to the left half-plane as the weights of the 
connections are varied. It turns out that r admits a simple interpretation in terms of the topology of the 
graph via the Mayer- Vietoris sequence. This, in turn, will provide a great deal of insight into the physics of 
the problem, allowing one to identify important structures in the network. We begin with a definition: 

Definition 2.17 Consider the following map dfrom i?i(T), the space of cycles in the graph, to R N . For each closed 
cycle 7 € -ffi(T) we associate the following vector v 1 — d(-y) e R N : 

• For each time time the cycle enters vertex i from a negatively weighted edge and exits through a positively 
weighted edge, v] increases by one. 

• For each time time the cycle enters vertex i from a positively weighted edge and exits through a negatively 
weighted edge, v] decreases by one. 

• For vertices not on the cycle, or for vertices where the cycle enters and exits through edges of like weights, v 1 is 
zero. 

It is clear that this map is linear, and that the image is an additive group. It should also be clear that 
cycles that remain entirely in edges of one type are in the kernel of this map and that only cycles with both 
types of edges give rise to nontrivial vectors v 1 . We refer to these cycles as "cycles of mixed type". 

Example 2.18 The following illustrates this map for two different graphs. The first graph has two cycles of mixed 
type. The first (71) is 3 — > 1 — > 2 — > 3 and the second (72) is 3 — > 5 — > 4 — > 3. This gives the additive group as linear 
combinations of v 11 = (—1,1,0, 0, 0) and v 12 = (0, 0, 0, 1,-1). The second graph has only one circuit of mixed type. 
The cycles 4^2^1^4 and 4->3->2^4 sum to 4 ^ 3 ^ 2 ^ 1 -> 4. Since this is a circuit of all positive 
edges it is in the kernel of d and thus the second circuit is the inverse of the first. This gives the additive group as 
vectors of the form k(0, 1, 0, —1). 
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Figure 2: The two graphs referenced in Example 



2.18 



Lemma 2.19 The flexibility 

T = iv + i-c(r + )-c(r_) 

is equal to the the dimension of the group of cycles of mixed type. 

Proof. This is a straightforward application of the Mayer- Vietoris sequence [34. §25]. We have the exact 
sequence 

Hi(r) -4 h q (t + n r_) A h (t + ) © ff (r_) A ff (r). 

The graph is assumed to be connected so dxm(H (T)) = 1, dim(iJ (r + ) © H Q (T_)) = c(T + ) + c(T_) and, 
since T + and T_ have no common edges, dim(i?o (T+ n T_ ) ) = iV. The exactness implies that dim(ker(/3) ) = 
dim(im(a)) = c(T + ) + c(T_) - 1 and dim(im(9)) = dim(ker(a)) = N - (c(T + ) + c(T_) - 1) = N + 1 - 
c(r_|_) — c(T_). Note that the last quantity is exactly the flexibility. Thus the flexibility is equal to the number 
of linearly independent cycles of mixed type. 

■ 

This construction shows that the dimension of the group of mixed cycles is equal to the number of 
eigenvalues which cross from the left to the right-half plane. It would be more satisfying to give an explicit 
bijection between the basis of mixed cycles in the graph and the eigenvalues which cross over, or the associ- 
ated eigenspaces. Of course we cannot really expect a bijection, since the eigenspace is an analytical object 
that depends sensitively on the edge weights, whereas the group of mixed cycles is topological and doesn't 
depend on the edge weights. Nevertheless, there is a sort of topological stand-in for the eigenspace which 
nicely characterizes the modes which do not cross over. First we make a definition. 

Definition 2.20 Let Shea subspace (chosen independently of the weights 7y) and P s be the orthogonal projection 
onto the subspace S. The subspace S is said to be a subspace of fixed index if the matrix 

P S C :S^S 

has the same index regardless of the choice of weights. A subspace of fixed index is maximal if 

dim(S) = c(r+) + c(T_) - 1. 

It is a simple consequence of the Courant minimax principle that the projection of a symmetric matrix 
onto a subspace cannot have more positive or more negative eigenvalues than the original operator. We 
know from theorem | |2.8) that for appropriate choices of the weights the matrix can have as few as c(r + ) — 1 
negative eigenvalues and (for a different choice of weights) c(T_) — 1 positive eigenvalues. It thus follows 
that a subspace of fixed index cannot have more than c(T+) — 1 negative eigenvalues, c(T_) — 1 positive 
eigenvalues and one zero eigenvalue, and that the maximum possible dimension of a fixed subspace is 
c(r + ) + c(T_) — 1. However it is not clear that one can actually have a maximal subspace of fixed index 
or, for that matter, any non-trivial subspace of fixed index at all. We conclude this section by showing that 
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there is always a maximal subspace of fixed index that has a natural topological construction. In essence 
this subspace gives a natural (orthogonal) decomposition into modes which do not have an eigenvalue 
crossing, and modes which do. 

First we need to define two complementary subspaces, one of which will be a maximal subspace of fixed 
index. The construction of these subspaces is similar in spirit to the construction of the cut-space and cycle- 
space from algebraic topology although the cut-space and cycle-space are subspaces of the vector space 
over the edge set, not the vector space over the vertices. For a nice description of the cut- and cycle-space 
and some applications to the theory of electrical networks see the paper of Bryant [51. 

Definition 2.21 Define S{ TCC to be the following subspace ofR N : Given a basis for the mixed cycles {7i}[ =1 let 
SWcc = span({v 7i }[ =1 .) 

Define S& xe d as follows: let T +i be the i th component ofT + , and let v l ' + be the characteristic vector ofT +i : 



i, 3 e r +! j, 
o, i0r +>i . 



Similarly, let v 1 ' be the characteristic vector o/F_ i . Then 

Lemma 2.22 The subspaces S{ ree and Sfi xe d are orthogonal complements: 

"Sfixcd = "Sfree 

Proof. First we check that dirn(Sg xe d) = c (r+) + c(r_) — 1. Note that the set of vectors {v 1 ^} is linearly 
independent, as is the set of vectors {v l '~}. However U {v 1 ,- } is not a linearly independent set, as 

one has 

c (r+) c (r_) 

E E = (1,1,1,.-., 1). 

i=l i=l 

We claim that this is the only relation. To see this note that we have the following identity for subspaces 
Q,R 

dim(span(Q U R)) = dim(<3) + dim(-R) - dim(span(Q n R)). 

Next note that span({??' + }) n span({??' _ }) consists of all vectors that are constant on components of r + and 
constant on components of T_ . Since T is connected this means that these vectors must be constant on all of 
r, and are thus proportional to (1, 1, 1, . . . , 1). It is clear that a basis for Sfree is given by any c(T+) + c(r_) — 1 
vectors from U {v 1 '^}. 

Next note that if 7 is a mixed cycle, with v 1 the corresponding vector, and v is constant on a component 
r +)i , then 

(if, if) = 0. 

To see this first note that the number of times the cycle 7 enters r + i must equal the number of times it 
leaves T + When it enters and leaves r +! ; it must do so through a negative edge. Each time it enters T + ; 
gives a +1 in some entry of v 1 , and each time it leaves gives a —1 entry. Thus (iP , v) is the sum of an equal 
number of +1 and —1 entries and is therefore zero. 

Since Sgxed an d Sfrec are orthogonal and have complementary dimensions they are orthogonal comple- 
ments of one another. ■ 

Definition 2.23 We define the following subspaces of ^fi xe d 

• Sfixed = G span^+l^ and (w, (1,1,..., 1)> = 0} 

. = {w\w e span^-}^ and (w, (1,1,..., 1)> = 0} 
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. Sg xed = span((l,l,l,...,l)) 

Lemma 2.24 The subspaces S^Z'° are C-orthogonal. Specifically, what we mean by this is if v and w are chosen 
from two different subspaces of S+ xcd , S nxcd , Sg xod , then 

(v,£w) = 0. 

Proof. S^ Kcd c ker(£), so the fact that any inner product involving that subspace is obvious. Thus we 
restrict our attention to S^ xcd . By definition v 1 '~ is constant on the component IV - and zero off of this 
component. Let dTy- denote the set of vertices which are not in IV but which are connected to it by a 
positive edge: the nearest neighbors of the component. By direct computation it is easy to see that £v l ~ 
takes the following form: 

( -Efeearv,. 7fc,j 3 G r ;',- 
\ £feerv_ 7k,j 3 € 9IV,_. 



(Cv 1 ■-) 



Thus Cv 1 '~ is a linear combination of vectors that are 1 on some vertex in IV and —1 on some vertex not 
in Ti> t - but connected to it by a positive edge. Each of these vectors is necessarily orthogonal to since 
v l,+ is constant on components of T+. ■ 
We need one more lemma to prove the final result, the well-known Sylvester theorem 

Lemma 2.25 (Sylvester's law of Inertia) If A is a square matrix and S a square non-singular matrix then A and 
B = AS have the same index. The matrices A and B are said to be Sylvester equivalent. 

Proof. This is an old result and many proofs are known. We include a short one here for complete- 
ness. Consider the one-parameter family of Hermitian matrices A(s) = (cos(s)J — isin(s)B')A(cos(s)I + 
i sin(s)ST). The matrices (cos(s)J— i sin(s)S^) and (cos(s)J-H sin(s)i?^) are both invertible, so dim(ker(A(s)) 
is independent of s, and there are no eigenvalue crossings. This gives a homotopy from A(Q) = A to 
A(^) = B^ AB that does not change the index, so the indices of the two must be equal. ■ 

Proposition 2.26 The subspace Sfixed is a maximal subspace of fixed index: The operator Ps tiyccd C (considered as an 
operator from S^cd to SWcd) has index 

md(p S!ixod c) = (c(r+) - 1, i, C (r_) - 1) 

independent of the choice of edge weights. 



Proof. Using Lemma 2.24 it follows that Ps liycDd C has the following block structure: 

P Sfixod £P Sflxod =£-®£o®£+ 

where C + = P Q + C, £_ = Pa- £ and £ = P q0 £ = 0. 

Consider the block £_ © £q. This arises by orthogonal projection of £ onto the set of vectors constant 
on r_, and thus this matrix is Sylvester equivalent to the following graph Laplacian: one contracts on the 
negative edges, giving a graph with vertices corresponding to the components of T_ . This is a standard 
graph Laplacian on a connected graph with c(T_ ) vertices and thus has c(T_ ) — 1 negative eigenvalues and 
one zero eigenvalue. Similarly £q © £ + is Sylvester equivalent to the negative of the graph Laplacian given 
by contracting on the positive edges. This has, by the same argument, c(T + ) — 1 positive eigenvalues and 
one zero eigenvalue. The zero eigenvalue is obviously counted twice in this argument, giving 

ind (Ps !M £P Srix J - (c(T_) - 1, 1, c(r+) - 1). 
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3 Refinements 



In this section, we present some refinements of Theorem 2.8 In Section 3.1 we present the Deletion- 
Contraction Theorem, which allows us to obtain recursive formulas for the A4(T(t)) polynomial in terms 
of "smaller" graphs obtained by deleting and contracting edges; this allows us to give a precise character- 



ization of the bifurcation structure of C(T(t)) in many cases. In Section 3.2 we discuss the asymptotics of 
the individual eigenvalues of C(T(t)) in the limits t — > 0, oo. Recall that in the statement and proof of The- 
orem 2.8 we discuss the number of eigenvalues in the left- and right-hand half-planes in these two limits; 



here we give more precise statements of the locations of these eigenvalues. 
3.1 Deletion-contraction theorem 

Definition 3.1 Let T = (V, E) be a weighted multigraph (loops and multiple edges allowed), and e € E(T) an edge. 
Let T/ e denote the graph obtained by removing edge e: the graph with vertex set V(T) and edge set E(T) — e.Ife is 
an edge which is not a loop let T. e denote the graph obtained by contracting on the edge, which is obtained as follows: 

• If edge e connects vertices v\ and v 2 then v% and v 2 are identified as a single vertex v*. 

• Each edge connecting a vertex to v\ or v 2 becomes an edge connecting that vertex to v* 

• Edge e is removed from the edge set. 

Note that ifT is a simple graph and v\ and v 2 are part of a triangle then the contracted graph T e will have multiple 
edges. Similarly if there are multiple edges connecting vi and i> 2 then F. e will have loops. 

We can now state the Deletion-Contraction Theorem. 

Theorem 3.2 (Deletion-Contraction theorem) Let The a weighted multigraph with e e E(T). Then M(T) can 
be computed by applying the following rules 

• If e is not a loop then M(T) = M(T/ e ) + 7 e ./Vf (T. e ). 

• Ife is a loop then M (T) = M(T /e ). 

• IfTis disconnected then M(T) =0 

Proof. The deletion-contraction recursion is well-known — see Chapter 13.2 of the text of Godsil and 
Royle 1 12 ] for one proof. ■ 
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Figure 3: The graphs T; T \ e; T.e, where e is chosen to be the edge with negative weight. 

Example 3.3 As an example, let us consider the graphs given in Figure^(here, think of the tasa symbol to separate 
out the terms containing the special edge). First, notice that V has three spanning trees not containing e, of weights 
8, 12, 24, and five spanning trees containing e, of weights — ICtt, — 15£, — 20t, — 30i, — 40£. By Lemma ^6\ 

M(T) =44- 115*. 
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Notice that every spanning tree ofT that does not contain e is also a spanning tree ofT \ e, and thus 

M(T\e) = 44. 

On the other hand, T.e has three spanning trees of weights 2, 7, 14. Therefore we have 

M(T.e) = 23. 

And we see that 

M(T) = 44 - 115* = 44 + (-5t)(23) = M(T /e ) - j e M(T. e ). 

In the context of dynamical systems, an important distinction to be made is that between graphs for 
which n + (C(T)) = and those for which n + (£(T)) > 0; the former are called stable and the latter unstable. 
This notion comes from the fact that if we consider the ordinary differential equation 

x = C(G)x, 

then the origin is stable to perturbations iff n + (£(G)) = 0; if not, then perturbations move away from the 
origin at an exponential rate. 

Moreover, as proved above, when we consider the homotopy T(t) = T + + tr_, the function n + (£(T(t))) 
is a non-decreasing function of t. From this, it follows that if we define 

t*(r) :=supn+(£(r(t))) = 0, 
t>o 

then for t < t* (T), the Laplacian is stable, and for t > t* (T), the Laplacian is unstable — in short, it undergoes 
a dynamical bifurcation. 

Two facts follow immediately from Theorem 2.8 first, that t*(T) > if and only if T + is connected, and 
t*(r) < oo if and only if T_ ^ 0. We will concentrate the most in what follows on the case of t* (T) 6 (0, oo). 

We also point out that rigid graphs undergo no bifurcation; in fact n + (T(t)) is constant on {t > 0}. From 
the bifurcation point of view, then, these are the least interesting cases. 

We now work out some special cases. 

3.1.1 One negative edge 

Let us first consider the case where T + is connected, and there is exactly one negative edge in T_, call it e. 
Then we have 

M(T(t))=M(T + )-t\ le \M(T.e). 

This means that 

■* = -M(r+) 

he\M(T.eY 

Notice that this formula gives t* > 0; since T + and T.e have non-negative entries, they will have opposite 
signs, since their dimensions differ by one (q.v. | |2.3[ l). 

Intuitively, we expect that the denser a graph is, the less powerful each individual edge would be. We 
present a couple of examples. 



Example 3.4 (Ring graph) Consider the case of the ring graph Rn, i.e. V(Rn) = {1, . . . , N} and 

(RN)ij 



1, j = i±l (mod N), 
0, else. 



We also define the path graph = Rn \ {1, N}. Choose any edge of Rn and flip its sign to minus one. According 
to the Proposition above, the homotopy C(T(t)) will lose stability at 

(k = M(T\e) _ M(P N ) 



M(T.e) M(R N -x) 
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It is easy to see that M(Rn) = N,M.(Pn) = 1/ so we have 

1 

t* = 



N — 1 



Example 3.5 (Complete graph) Consider the complete graph Kjq. Again choose any edge, then the deletion opera- 
tors give graphs OGn,DRn~i- OGn is the complete graph minus one edge; DR^ is a complete graph on N nodes 
with all of the edges going to one distinguished vertex having twice the weight. 
Using standard counting arguments, we see that 

M{OG N ) = (N- 2)N N - 3 , M{DR N ) = 2{N + 1) N ~ 2 , 

giving 

M{OG N ) _ (N - 2)N N ~ 3 _ N - 2 



M(DR N ^) ' 2N N - 3 2 ' 

Of course, this can be obtained by other means; in fact, one can compute that the eigenvalues ofT(t) are 

{0} U {-iV}^- 2 ) U {2t + (N - 2)}, 

and the computation would also follow from this. 

Conjecture 3.6 These are the extreme cases; for any graph T with T + connected and \E(T_)\ = 1, we have 

t*(T)e[(N-l)-\(N-2)/2}. 
Proposition 3.7 Assume that |i?(r_)| = 2 and that the two edges e, / e E(T-) do not share a vertex. Then 
M(T(t)) = t 2 | 7e7/ | M(T. eJ ) - t (| 7e | M(T. e \ f ) + | 7/ | M(T v . f )) +M(T VXf ), 



where we have defined all of these terms in Definition 3.1 Also, t* (T) is the minimal positive root of this polynomial. 



Proof. Let r be any signed graph; from Theorem 3.2 we have 

M(T)= le M(r. e )+M(T /e ), 

and 

M(T. e ) = lf M((T. e ) f ) +M((T. e ) /f ), M(T/ e ) = lf M((T /e ) f ) + M((T /e ) /f ), 
and iterating \3.2\ gives 

M(T) = 7e7/r. e ./ + 7eM(r. Ae ) + 7 f M(T, e \ f ) + M(T V \ f ). 
If the edges e and / are in E(T-), then 7e = — t | 7e | , 7 / = —t | 7 /|, and the result follows. 



Remark 3.8 Note that this polynomial has positive roots by Theorem 2.8 Consider, for example, an unweighted 
graph; this implies (q.v. Remark 2.14\ that 

(M(T.e \ /) + M(T \ e.f)f > 4M(T.e.f)M(T \ e \ /). 
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3.2 Detailed Eigenvalue Asymptotics 



We now consider more detailed asymptotics of the eigenvalue spectrum in the limits in which the strength 
of the negative edges is much weaker or much stronger than the strength of the positive edges. Specifically 
we consider the one-parameter family of graph Laplacians 

in the limits t -> 0+ and t -> oo. The spectrum splits naturally into two parts, which correspond to the 
eigenvalues of graph Laplacians on the deleted and contracted graphs. This is an analog on the level of 
the spectrum of the contraction-deletion algorithm for computing the crossing polynomial: the crossing 
polynomial is given by the sum of the crossing polynomials for the contracted and deleted graphs, while 
the spectrum is given (asymptotically!) by the union of the deleted and contracted graphs. 

Theorem 3.9 Suppose that t is large and positive. Then C(t) has exactly N — c(T_ ) negative eigenvalues, c(r_ ) — 1 
positive eigenvalues and one zero eigenvalue. If we take the convention that eigenvalues are numbered in decreasing 
order then to leading order in t the c(T_ ) — 1 positive eigenvalues are given by 

A,(£ r ) = t\ l (£ r _ ) + 0(1) i e {1 . . . c(T_) - 1}, 

The N - c(T_) negative eigenvalues are given to leading order by 

\(Cr) = Ai(r r .-) + o(l) i € {N . . . c(T_) + 1}, 

where (Lr— is the graph formed by contracting on the negative edges. Here Xj (£r— ) are solutions to 

Ly-v = XSv 

where L T - is the graph Laplacian formed by contracting on the negative edges and S the contracted inner product: 
the diagonal matrix with entries S li = |"V"(r ,»)|. 

Proof. The proof follows in a straightforward way from perturbation theory for eigenvalues of a 
symmetric matrix. The full graph Laplacian can be written 

so it suffices to understand the eigenvalues of £(T_) + t~ 1 £(T + ) for t large. The spectrum of £(T_) consists 
of c(r_ ) zero eigenvalues and N— c(r_ ) positive eigenvalues. For the non-zero eigenvalues straightforward 
eigenvalue perturbation theory gives the asymptotic above. 

To understand how the c(r_)-dimensional kernel breaks under perturbation we must do a degener- 
ate perturbation theory calculation. Well-known results (again see Kato) show that to leading order the 
eigenvalues are given by the eigenvalues of the reduced matrix 

^kor(£_)£r + -Pkcr(£_) 

where Pkcr(£_) is the orthogonal projection onto the kernel of £_. It is straightforward to compute ker(£_): 
it consists of vectors that are constant on components of We can identify a vector w e kcr(£_) with a 
vector w € M c ( r ~' by the following rule: if Wi = a for all vertices in component r_, j then Wj = a. Under 
this identification the natural inner product on R N maps to the inner product 

C (r_) 

(v,w)= J2 

i=\ 

In other words the there is one entry per component of T_, and the inner product is diagonal with weights 
given by the number of vertices in the corresponding component of IV It is straightforward to see that the 
matrix 

^ker(£_)-Cr + -Pker(£_) 

is exactly the Laplace matrix obtained by contracting on the negative edges of the graph, completing the 
proof. ■ 
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Definition 3.10 Given a matrix A, we define Spec* (A) as the eigenvalues of the matrix A restricted to the subspace 
(1, 1,1,..., I) -1 of mean zero vectors. 

Remark 3.11 The previous theorem can be written in the following compact way 



Spec*(£(t)) 



Spec*(r_)USpec*(r._) 
Spec*(r+) U Spec*(r.+) 



+oo 
0+ 



with the understanding that the eigenvalues for the contracted graph are taken with respect to the natural inner prod- 
uct S. This shows that there is an (approximate) contraction-deletion relation at the level of the spectrum analogous 
to the contraction-deletion relation satisfied by the crossing polynomial. 

Example 3.12 We consider the following graph, where all positive (solid) edges are weighted +1 and all negative 
(dashed) edges weighted — 1; The full graph has five positive edges, seven negative edges and nine vertices. The 




A 



o 



D 



5 

o™ 




Figure 4: The graph T; the graph T_; the graph T._ 

subgraph T_ consisting of only the negative links has three components. Component A consists of vertices 1,2 and 3 
and connecting edges , component B consists of vertices 4 and 5 and the connecting edge, and Component C consists 
of vertices 6, 7, 8, and 9 and connecting edges. There are N — c(r_) =9 — 3 = 6 non-zero eigenvalues corresponding 
to the graph Laplacian associated with the negative edges. The non-zero eigenvalues associated to component A are 1 
and 3, to component B is 2 and to component C are 4, 2 and 2. This gives six eigenvalues that grow linearly: 



A 9 
A 8 
Ay 
A 6 
A 5 
A 4 



4t + 0(l) 
3t + 0(l) 
2t + 0(l) 
2< + 0(l) 
2t + 0(l) 
lt + Ofl) 



If one contracts on all of the dashed edges the full graph reduces to the three cycle, with one vertex corresponding to 
each component ofT— There are two edges between components A and B, two between components B and C, and one 
between A and C, so these edges are weighted accordingly. The norm is contracted as well, and the norm over the new 
vertex space can be written as 

\\v\\ 2 =i7*5w 

where S is the matrix 

/ 3 \ 
S= 2 
\ 4 / 
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The diagonal entries reflect the fact that the components have three, two and four vertices respectively. Thus the 
eigenvalue problem becomes 



giving the negative eigenvalues as 
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1 



33- 151 +0(l/t) 
-V33- 15) +0(l/t) 



Aa(<) 
Ai(t)« 

Finally the flexibility is equal to t(T) = 10 — 3 — 4 = 3, so there are three eigenvalue crossings. It is straightforward 
though tedious to compute that the crossing polynomial is given by Pr(t) = 171t 3 — 702i 4 + 828i 5 — 288i 6 . The 
non-zero roots occur att~ .43, t ~ .90, t ~ 1.55. 

Some numerical results are shown in Figure . The first plot shows a graph of Xi(t)/tfor i = 3 . . . 9 and t e (0, 8). 
It seems clear that the scaled eigenvalues are converging to the correct values. The second plot shows a plot of the 
(unsealed) negative eigenvalues for t £ (0,8). Again it is clear that they are converging to — 4(\/33 + 15) and 
— I (15 — a/33) respectively. One can also see that there are three eigenvalue crossings at the correct t values. 




Figure 5: Eigenvalues of T(t) as a function of t, where T(t) is defined in Example 3.12 



3.3 Comparison to the Gershgorin Theorem 

In applied mathematics and numerical analysis it is often necessary to estimate the locations of the eigen- 
values of a matrix or linear operator. One simple and very widely used tool for this is the Gershgorin disc 
theorem, which says the following: 



Theorem 3.13 (Gershgorin) Given a matrix M with entries My define the following n closed disks: 

Di = {z\\z - Ma 



<£lMii>. 



Then the eigenvalues of M lie in the union of the disks, spec(M) € U" =1 -Dj. 

The fact that the disks are defined using the edge weights might make one think that the Gershgorin 
theorem might give more information on the signs of the eigenvalues than the purely topological argu- 
ments, but this is not, in fact, the case. It is not difficult to see that, for a non-trivial signed Laplacian the 



Gershgorin theorem always gives results that are strictly worse than those given by Theorem 2.8 

First note that if the matrix M represents a graph Laplacian then there are three kinds of discs. If all 
of the edges emanating from the vertex are all positive (resp. all negative) then the associated disc lies 
in the closed left half -plane (resp. right half -plane) and is tangent to the origin and the corresponding 
eigenvalue is non-positive (resp. non-negative). If the vertex has edges of both signs then the origin lies in 
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the interior of the disk and the sign of the eigenvalue cannot be determined. It is these discs that correspond 
to eigenvalues whose sign cannot be determined. The main observation in this section is that the number 
of such discs is always strictly larger than r. 

Proposition 3.14 Suppose the graph T is connected and contains edges of both signs. Let n be the number of Gersh- 
gorin discs for which the origin lies in the interior. Then n > r + 1. 

Proof. This follows immediately from the topological characterization of the number of of modes for which 
there is an eigenvalue crossing. The results of the previous section show that these can be associated with 
mixed cycles. By the construction given in that section we associate to each of these cycles a vector that has 
+1 for each time the cycle enters a vertex from a positive edge and leaves it by a negative edge, and —1 for 
each time the cycle enters a vertex from a negative edge and leaves it by a positive edge. These vectors can 
only have non-zero entries in vertices which have both types of edge, so they obviously lie in a subspace 
isomorphic to R". These vectors are necessarily orthogonal to (1, 1, 1, . . . , 1), so there can be at most n — 1 
linearly independent ones. ■ 

Remark 3.15 Note that n = r + 1 can be achieved — one example is when c(T+) = 1 and c(T_) = 1, when 
t = N — 1 and n = N. It can also happen that n is much larger than r. Consider, for example, an even cycle with 
edges of alternating sign. In this case n = N — all vertices have edges of both type — while r = 1 since there is only 
one linearly independent loop. 

Example 3.16 The graph depicted in Figure^has a flexibility ofr = 9 + 1 — 4 — 3 = 3. This graph has eight vertices 
which have edges of both signs, and thus eight Gershgorin discs that contain the origin in the interior. 

In the second graph in figure [I] there are three vertices that have edges of both types, and thus three Gershgorin 
discs that contain the origin as an interior point. The flexibility of this graph, however, is zero so that the number of 
positive, negative and zero eigenvalues is fixed and does not vary with the edge weights. 

4 Applications and numerical computations 

4.1 Random graphs and bifurcations 

It was shown in Section |3T| that for any signed T, 

t* := sup{t : n+(r(t)) = 0} 
t>o 

satisfies t* G [0, oo], and that t* = if and only if T + is not connected, and t* = oo if and only if T_ = 0. 
Thus we can think of t* as a map from the set of finite graphs to [0, oo]. For any collection Q of graphs and 
a probability measure P on Q, this induces a random variable T: Q — > [0, oo]. 

We present some numerically-computed distributions in Figure[6]and[7]below. The random ensemble of 
matrices is a signed generalization of the classical Erdos-Renyi random graph G(n,p) whose specific defi- 
nition is given as follows: we fix N and p± € [0, 1]. For each 1 < i < j < N, choose Xij, Yij independently, 
with 

F(Xij = l)= P+ , P(X y = 0) = l-p+, P(i; J = l)=p_, P(Yy = 0) = l-p_, 
and we set 7^ = — . From this, we have 

P(7ij = +1) = P+(l - P-), P(7« = -1) = P-(l -P+), 

and jij — otherwise, and of course the 7^ are independent as well. We then set 7^ with i > j by 
symmetry. This gives a random distribution on the set of symmetric graphs with N vertices. We condition 
on r + being connected and T_ 7^ 0, and then compute the distribution of t* over this ensemble. 

We performed a series of numerical experiments on these random variables, both for N = 10 ("small 
matrices") and N — 50 ("large matrices"), the numbers 10 and 50 being chosen arbitrarily. In each of these 
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N = 1 0, p_ = 2.00e-01 N = 10, p = 2.00e-01 




Figure 6: Each of the rows corresponds to a fixed N; the top row is N = 10 and the bottom row is N = 50. In 
each row, in the left panel, we plot in solid blue the mean and standard deviation for ensembles of random 
matrices where we vary p + , and in light blue, we plot the value of t* for each individual matrix. In the right 
panel, we have plotted the distributions for all values of p+, but rescaled to have mean zero and variance 
one. For the case of N = 50, we compare these rescaled distributions to the standard Gaussian. 



cases, we fixed p_ for all simulations (we chose p_ = 0.20 for TV = 10 and p_ = 0.10 for N = 50) and 
varied p+ over a range. For each choice of p±, we chose a random ensemble of 10 4 matrices using the rules 
above, and computed t* for each of the matrices in the ensemble using a bisection method. We present the 
findings in Figures [6] and [7] We first observe that t* tends to increase as p+ increases, which makes sense: 
adding more positive edges will make the Laplacian more stable, and higher p + values tend to give more 
positive edges. What is perhaps surprising is that the ensemble mean is very close to a linear function of 
p+ (at least for a certain range for N = 10, and for all p + for N = 50). Moreover, we find that there seems 
to be a universal scaling distribution for the different values of p+ ; in each case, what we do is consider the 
distribution of t* for each choice of p±, then normalize this distribution to have mean zero and variance 
one, and plot these on top of each other. For the small (N = 10) case, the rescaled distributions overlap well 
and resemble a lognormal plot; for the N = 50 case, the rescaled distributions overlap well and resemble a 
normal plot. 

To make these observations more precise, in each case we chose p+ = 0.45 and made a QQ-plot of the 
ensemble distribution versus either a lognormal (for N = 10) or a normal (for N = 50), and these match 
quite well. See Figure [8] 



4.2 Feuds in social networks 

We consider two datasets from social networks and compute the flexibility of the graphs and, in one case, 
the bifurcations. Our computations were facilitated by the matlab_bgl library^] 

The first dataset we consider is from Read |37|, and represents sympathetic and antagonistic relation- 
ships amongst sixteen sub-tribes of the Gahuku-Gama people in the highlands of New Guinea: Gaveve, 
Kotuni,Gama, Nagamidzhuha, Seu've, Kohika, Notohana, Uheto, Nagamiza, Masilakidzuha, Asarodzuha, 
Gahuku, Gehamo, Ove, Ukudzuha and Alikadzuha. Warfare between subtribes in this society was very 
common. While temporary alliances are common this dataset represents traditional relationships between 
subtribes. The positive, or hina, edges represent subtribes that are traditionally close politically. Warfare be- 

' http : / / www . mathworks . com/mat labcentral/f ileexchange/10 922 
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N = 50, p_ = 1 .00e-01 N = 50, p_ = 1 .00e-01 




Figure 7: In the left panel, in solid blue we plot the mean and standard deviation for an ensemble of random 
matrices with a given p + , p_ = 0.1, and N = 50, versus p + . In light blue, we plot the value of t* for each 
individual matrix. In the right, we have plotted the distributions for all values of p+, but rescaled to have 
mean zero and variance one. We compare with the Gaussian in red. 

tween these subtribes occurs but is limited. It is understood to be short term condition and is often resolved 
by payment of blood money or other concessions. The negative, or rova, edges represent relations between 
subtribes that are traditionally antagonistic. Warfare between these subtribes is much less constrained. This 
has become a somewhat popular data set to analyze: see the pioneering work of Hage and Harary 1 16] and 
the recent work of Kunegis et. al[28|. The graph has c(T + ) = 2andc(T_) = 3,givingr = 12. The first com- 
ponent (A) of r + consists of a collection of four tribes (Gaveve, Kotuni,Gama, Nagamidzhuha; here num- 
bered 1-4) all of whom have friendly relations. The second component (B) consists of the remaining twelve 
tribes (Seu've, Kohika, Notohana, Uheto, Nagamiza, Masilakidzuha, Asarodzuha, Gahuku, Gehamo, Ove, 
Ukudzuha and Alikadzuha; here numbered 5-16), all of whom are connected by at least one chain of sympa- 
thetic relationships. Relations are more complicated within component (B) than in component A but there 
are two main features to be noted. The first is that removing the Masilakidzuha (tribe number 10) splits 
this component into two subcomponents, the first subcomponent (Bl) consisting of the Seu've, Kohika, 
Notohana, Uheto, Nagamiza and the second subcomponent (B2) consisting of the Asarodzuha, Gahuku, 
Gehamo, Ove, Ukudzuha and Alikadzuha. There are numerous antagonistic relationships between tribes 
in subcomponents Bl and subcomponent B2 but there are no antagonistic relationships within subcompo- 
nents Bl or B2. This suggests the possibility of a rift developing within component B. 

The second dataset is from the Slashdot Zoo. Slashdot[jis a user-run website where links are submitted 
and voted on by the userbase, and, furthermore, users can make comments on the links and these com- 
ments are also voted on by the individual users. A significant amount of discussion occurs on this website, 
sometimes friendly and sometimes fractious, and thus it is not hard to imagine that connections, both pos- 
itive and negative, form between users. Each user is allowed to tag other users in the database as a "fan" 
or a "foe", i.e. if user i likes the types of comments made by user j, user i has the possibility to tag user j 
and become a "fan". If this occurs, then in the Slashdot Zoo [^database, the edge (i, j) is given weight +1. 
Similarly, user i can be a "foe" of user j and this adds a —1 on edge This dataset contained connection 
data on 82,144 users, with 549,202 edges in the graph. This network is not a priori symmetric, since the 
fan/foe operations have directionality, so we imposed symmetry. Given users i and j, if i and j are both 

' http : //www . slashdot ■ com| 

f http : / /konect . uni-koblenz . de /networks / s lashdot - zool dump of userbase, May 2009 
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QQ Plot, N=1 versus longnormal QQ Plot, N=50 versus Normal 
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Figure 8: QQ-plots of ensemble distributions. In the left frame, we are plotting a QQ-plot of the ensemble 
distribution N = 10, p- = 0.20, p+ = 0.45 versus a lognormal; in the right frame, we plot the ensemble 
distribution for TV = 50, p_ = 0.10, p + = 0.45 versus a normal. 

fans of each other, or i is a fan to j and j is neutral to i, then we placed a +1 in edge (i, j). Similarly for foes, 
if i and j are both foes, or i is a foe to j and j is neutral to i, then we placed a —1 in edge In short, 

we allowed simply extended unidirectional relationships to be bidirectional as long as the other direction 
was neutral. The only case to think about is what one should do if the two directions are of opposite sign, 
i.e. if i was a fan of j and j a foe of i; in this case we decided to assume that the relationship canceled and 
placed a on edge As one can imagine, this is relatively rare, and this happened only 1,949 times in 

this dataset. 

Thus, in both cases, we are considering a network with positive and negative weights. Since it has a 
smaller vertex set, the PNG dataset was easier to analyze. We present this data in Figure[9] Since c(T + ) = 2, 
and t = 12 the results of Theorem |Z8] imply that we will have we will have 

lim n + (T(t)) = 1, lim n+(T(t)) = 13. 

t-tQ+ t-¥oo 

A symbolic computation using Mathematica gives the crossing polynomial as 

M(t) = -45432223i 13 + 657635624* 12 - 4187415940t n + 15505043366i 10 - 37159886129t 9 

+ 60647687776i 8 - 6896052657K 7 + 54844706645i 6 - 3010376212U 5 + 11015925656* 4 
- 2508107376* 3 + 308319872t 2 - 14192640i. 

The twelve positive roots of the crossing polynomial range from t w .1 to t w 3. For small positive t we 
have one positive eigenvalue. This, of course, represents the mutual antagonism between the subtribes 
in component A and the subtribes in component B. The first bifurcation occurs around t w .1 when a 
secondary instability develops. Intuitively one expects that this corresponds to a split in component B, 
with antagonism between sub-components Bl and B2. This is confirmed by the numerics. If we compute 
the spectrum of the Laplacian at the first bifurcation point t w .988 there is a single positive eigenvalue with 
eigenvector 

v = (-.43, -.43, -.44, -.43, .14, .16, .16, .15, .16, .13, .16, .15, .14, .14, .13, .12)*. 

This obviously represents the mutual aggression between components A and B. There is a second linearly 
independent vector in the kernel of the Laplacian representing the emerging instability. This eigenvector is 
given by 

v = (.03, -.002, .01, -.003, .34, .42, .38, .30, .22, -.12, -.27, -.29, -.29, -.24, -.24, -.25)*. 
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We see that for this mode component A is only weakly involved - component A has only about 3% of 
the mass of the eigenvector - and the eigenvector clearly describes a rift in component B. The tribes in 
subcomponent Bl (subtribes 5-9: Kohika, Notohana, Seu've, Uheto and Nagamiza) move in one direction 
and the tribes in subcomponent B2 (subtribes 11-16: Asarodzuha, Ove, Gahuku, Gehamo, Ukudzuha and 
Alikadzuha) move in the opposing direction. The Masilakidzuha (subtribe 10) have somewhat stronger 
ties to subcomponent B2 than to subcomponent Bl (five sympathetic relationships with tribes of B2 vs. two 
sympathetic relationships with tribes in Bl) and so move with subcomponent B2. 

The next bifurcation, which occurs around t « .57, describes a somewhat less obvious conflict. The 
associated null-vector at the bifurcation point is 

v = (0.44, 0.29, -0.26, -0.37, -0.05, 0.16, -0.002, 0.19, -0.23, 0.024, 0.25, -0.04, 0.32, -0.29, -0.022, -0.40)* 

which describes a conflict with subtribes 1,2,6,8,11 and 13 (Gaveve, Kotuni, Notohana, Uheto, Asaro- 
dzuha and Gahuku) forming one faction and subtribes 3,4,9,14 and 16 (Gama, Nagamidzhuha, Nagamiza, 
Gehamo and Alikadzuha) forming another, with much smaller involvement of the remaining tribes. As 
t increases we see additional eigenvalue crossings, leading to additional conflicts, up to the maximum of 
thirteen unstable modes for t ps 3. Moreover, if we assume that friend and foe links are of equal strength, 
then we have n + (T) = n + (T(l)) = 7. 



a) 



§■ 8 

o 

i?io 




5 10 15 0.5 1 1.5 2 2.5 3 3.5 

tribal group number 1 



Figure 9: The PNG dataset from \37\. In frame (a), we are plotting the interaction matrix amongst the 16 
tribes. Blue pixels correspond to entries in T + , red to entries in T_. In frame (b), we plot n + (T(t)) as a 
function of t; we know that this is positive for all t > 0, since c(T + ) > 0, and we see that it increases to its 
maximum value around t = 3.0. 

For the Slashdot Zoo data set, we concentrated on the component of the Slashdot database that is 
friendly to CmdrTaco, the founder of the site and user number 1. More specifically, we considered only 
those users for which there existed a "friendly path" from that user to CmdrTaco. This subset was, not 
surprisingly, the largest connected component of the full network, and contains 23,514 users. We then con- 
sidered the subgraph of these 23,514 users to itself, and this is what we define V to be. This graph contains 
415,118 positive edges and 117,024 negative edges, and we computed components. By definition, c(T + ) = 1, 
and we compute that c(T_) = 10,327, giving T a flexibility index of 13,187. This is, admittedly, a one-off 
calculation, but we conjecture that data from social networks will show this pattern, that even amongst a 
group of "friends", or common "fans" of a particular user, there will be a large number of instabilities in 
exactly this manner. 
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